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Abstract It is shown that the initial singularities in spatially compact spacetimes with 

. spherical, plane or hyperbolic symmetry admitting a compact constant mean curvature 

Q>^ [ hypersurface are crushing singularities when the matter content of spacetime is described 

"q I by the Vlasov equation (coUisionless matter) or the wave equation (massless scalar field). 

O^' In the spherically symmetric case it is further shown that if the spacetime admits a maximal 

Vh . slice then there are crushing singularities both in the past and in the future. The essential 

properties of the matter models chosen are that their energy-momentum tensors satisfy 

certain inequalities and that they do not develop singularities in a given regular background 

/S ' spacetime. 



1. Introduction 

The nature of singularities in general solutions of Einstein's equations is still a matter 
about which very little is known. The best information on this which has been obtained 
up to now concerns the singularities in special classes of solutions defined by symmetry 
assumptions. The hope is that the insights obtained in solving these restricted problems 
will allow the symmetry assumptions to be progressively relaxed and so the study of 
singularities in solutions of the Einstein equations with various symmetries can be seen 
as a systematic approach to the general problem of understanding spacetime singularities. 
The following is intended as a contribution in this direction. 



* Present address 



At the moment the study of singularities in solutions of the vacuum Einstein equations 
is significantly further advanced than in the case of the Einstein equations coupled to 
matter. This paper is concerned mainly with non-vacuum spacetimes and helps to redress 
the balance a little. The spacetimes treated have spherical, plane or hyperbolic symmetry 
and a compact Cauchy hypersurface. (The reader wishing to gain some intuition for this 
class of spacetimes is referred to Appendix 2, where the vacuum solutions are determined 
explicitly.) Some of the results are for general matter models which are only restricted 
by some inequalities on the components of the energy-momentum tensor. However the 
main results require a matter model which is well-behaved in a certain sense. Roughly 
speaking, the matter fields should not develop singularities in a given regular spacetime. 
These results are worked out for two particular matter models, namely coUisionless matter 
and the massless scalar field (Theorems 3.2 and 3.3 respectively). 

There are some differences between the results for the different symmetry classes. 
For plane symmetry it is shown that, if the spacetime is the maximal globally hyperbolic 
development of constant mean curvature (CMC) initial data on a compact spacelike hyper- 
surface, there exists a foliation of compact hypersurfaces of constant mean curvature which 
covers either the past or the future of the initial hypersurface. If the mean curvature of the 
initial hypersurface is Hq then the range of the mean curvature of this foliation contains the 
interval (— oo, Hq] or [Hq, oo) respectively. The analogous result in the case of hyperbohc 
symmetry is proved under the additional assumption that the mass function, defined by 
equation (2.15), is positive on the initial hypersurface. In the spherically symmetric case 
the results are as follows. If the initial hypersurface is maximal (i.e. its mean curvature is 
zero) then it is shown that the spacetime can be covered by a CMC foliation where this 
time the range of the mean curvature is the whole real line. This should probably be true 
whatever the mean curvature of the initial hypersurface but all that could be shown is 
that there is a foliation by compact CMC hypersurfaces whose range includes all positive 
real numbers or all negative real numbers and this foliation will in general not cover the 
whole spacetime. It does however cover the part of the spacetime either to the past or 
to the future of the initial hypersurface. All these results are dependent on the choice of 
a well-behaved matter model, as indicated above. In the spherically symmetric case they 
are also dependent on the assumption that the topology of the Cauchy hypersurface is 
5"^ X 5"^ so that there is no centre of symmetry. The differences between the various cases 
are related to the fact that while an initially expanding solution with plane or hyperbolic 
symmetry can be expected to expand for ever a spherically symmetric solution can be 
expected to recoUapse. 

In [6] and [7] conjectures were formulated concerning the existence of CMC foliations. 
The results just discussed are closely related to conjecture 2.3 of [7] and conjecture C2 
of [6] for the class of spacetimes considered here. They show in particular that any of 
these spacetimes has a crushing singularity in at least one time direction (and in both 
time directions if it contains a maximal hypersurface). Recall that a crushing singularity 
is one where there exists a foliation on a neighbourhood of the singularity whose mean 
curvature tends uniformly to infinity as the singularity is approached [7]. A difficulty in 
studying spacetime singularities which may exist at the boundary of the maximal Cauchy 
development is that the definition of the maximal Cauchy development is so abstract. 



(Zorn's Lemma is used to show its existence.) Having a geometrically defined global time 
coordinate helps to make it more concrete. It allows the strong cosmic censorship hypoth- 
esis to be reformulated as a question on the global existence and asymptotic behaviour of 
solutions of a system of partial differential equations[6]. In this sense proving the existence 
of global CMC foliations in a given class of spacetimes is a first step towards proving strong 
cosmic censorship in this class. It also guarantees that a numerical calculation done using 
a CMC slicing is in principle capable of covering the whole spacetime. (The possibility of 
using a CMC slicing for numerical studies of cosmological spacetimes has been discussed 
in [9].) 

The most general results on the existence of CMC foliations in any class of spacetimes 
with compact Cauchy hypersurfaces are due to Isenberg and Moncrief[10]. They showed 
that Gowdy spacetimes (which are vacuum spacetimes with U{1) x U{1) symmetry) on 
the torus can be foliated globally by constant mean curvature slices. For spacetimes with 
matter less is known. In [13] it was shown that spatially homogeneous spacetimes with 
matter described by the Vlasov equation can be foliated by CMC hypersurfaces with the 
mean curvature ranging either over the whole interval (—00, 0) or the whole real line, 
depending on the symmetry type. In [15] a similar result was obtained for perfect fiuids 
and for non-interacting mixtures of perfect fiuids. 

There is another motivation for the work reported in this paper which has not yet 
been mentioned. When investigating the global structure of spacetimes it seems essential 
to introduce some auxiliary elements such as coordinates in order to 'find one's way in 
spacetime'. If the spacetimes being studied have symmetry properties then it may be 
advantageous to use elements which exist due to the symmetry. The problem with this, 
from the point of view of the general programme outlined above, is that this approach 
is less likely to produce techniques which can be used beyond the given symmetry type. 
The definition of a CMC hypersurface does not depend on any symmetry assumptions and 
there is at least a chance that global CMC foliations will exist in rather general spacetimes. 
One aim of this work was to see whether it was possible to prove something about the 
properties of singularities working with a CMC foliation from the beginning. This turned 
out to be the case. 

It is plausible that using coordinates adapted to a special symmetry type could be 
more efficient for proving sharper results than CMC slicing. This is confirmed by the 
fact that more precise information about the nature of the singularities in the spacetimes 
considered in this paper has been obtained in certain cases using other coordinate systems. 
In [14] the case of a plane symmetric scalar field was studied. It was shown that the initial 
singularity is a curvature singularity and a velocity-dominated singularity in that case. 
These results provide a model for what one would like to do more generally. For the case 
of collisionless matter Rein[12] has shown that under certain assumptions on the initial 
data the initial singularity is a curvature singularity and possesses some attributes of a 
velocity-dominated singularity. 

2. Analysis of the field equations 

The spacetimes of interest here are defined on manifolds of the form M = Hx S^ x F, 
where F is a compact orientable surface. Let p denote the projection of the universal 



cover F onto F. Let Qccp be a globally hyperbolic metric on M for which each submanifold 
{t} X S^ X F is a, Cauchy hypersurface. The spacetime (M, g^p) is called spherically 
symmetric if F = S'^ and the transformations of M induced by the standard action of 
SO{3) on S'^ are isometries which leave any matter fields invariant. In the case F = T^ the 
projection p induces a projection p of M = R x 5"^ x R^ onto M. Let (jctp be the pull-back 
of QafB by p. The spacetime is called plane symmetric if F = T^ and if the tranformations 
of M induced by the standard action of the Euclidean group E2 on R^ are isometries which 
leave (the pull-backs of) any matter fields invariant. In the case where the genus of F is 
greater than one F can be identified with the hyperbolic plane H^. A projection p can be 
defined as in the plane symmetric case. The spacetime is said to have hyperbolic symmetry 
if the genus of F is greater than one and the transformations of M = R x 5"^ x H^ induced 
by the action of the connected component of the identity of the isometry group of the 
hyperbolic plane are isometries which leave any matter fields invariant. In the spherically 
symmetric case SO{3) acts on 5"^ x 5"^ without fixed points. In other words there are no 
centres. The group SO (3) also acts on 5"^ with fixed points and this leads to a different class 
of spherically symmetric spacetimes. This other class will not be considered in this paper 
and so for clarity the spherically symmetric spacetimes considered here will be referred 
to as spherically symmetric spacetimes without centre. It will be convenient to refer to 
spherically symmetric spacetimes without centre and spacetimes with plane and hyperbolic 
symmetry collectively as surface symmetric spacetimes. The surfaces diffeomorphic to F 
which are defined by the product decomposition will be referred to as surfaces of symmetry. 
Now the Einstein equations for surface symmetric spacetimes will be analysed in a 
certain coordinate system. In order to ensure the existence of such a coordinate system 
it will be assumed that the spacetime possesses a Cauchy hypersurface of constant mean 
curvature which is symmetric, in the sense that it is a union of surfaces of symmetry. 
Lemma 2.1 Let {M,g) be a non-fiat surface symmetric spacetime having a symmetric 
constant mean curvature Cauchy hypersurface and satisfying the dominant energy and 
non-negative pressures conditions. Then in a neighbourhood of any point there exist local 
coordinates adapted to the product decomposition R x 5"^ x F such that the metric takes 
the form 

-a'^dt^ + A'^[{dx + /3dtf + a^rfE^]. (2.1) 

The functions a, (3 and A depend on t and x, a depends only on t and the metric dT,^ 
has constant curvature e. It may be assumed without loss of generality that e = 1 for 
spherical symmetry, e = for plane symmetry and e = — 1 for hyperbolic symmetry. The 
time coordinate t may be chosen so that the hypersurface where t has a given value has 
constant mean curvature equal to that value. 
Proof See Appendix 1. 

The neighbourhood in this lemma may be chosen to be the product IiX I2X U, where 
/i and I2 are intervals and U is an open subset of F. The interval /i will be denoted by 
(ti, ^2) in the following. By a slight extension of the usual notion of coordinates I2 can be 
taken to be the closed interval [0, 2tc], where it is understood that x = and a; = 27r are 
to be identified. The functions a, /3 and A are assumed to be functions which have C°° 
extensions to functions which are 27r-periodic in x. The coordinates can be chosen so that 
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/q /3(t, x)dx = for each t. The second fundamental form of each hypersurface t = const, 
can be written as 

A^{Kdx^-^{K-t)dJ:^) (2.2) 

for a function K{t, x) which has the same regularity properties as those demanded of ct, (3 
and A above. The Einstein equations will now be written out for the metric (2.1), making 
use of the variable K defined by (2.2). 

(Ai/2)" = _i a5/2[|(^ _ 1^)2 _ 2^2 ^ ;L67rp] + ieA-V2a-2 (2.3) 

a" + A-^A'a' = aA%{K - \tf + \e + Aii{p + trS)] - A^ (2.4) 

K' + ?,A-^A'K - A-^A't = SttjA (2.5) 

^' = -a-^dta + \a{?,K-t) (2.6) 

dta = a[-l3' + \a{?,K-t)\ (2.7) 

dtA = -aKA + (PA)' (2.8) 
dtK = (3K' - A-^a" + A'^A'a' 

+ a[-2A-^A" + 2A-^A'^ + Kt - 8nSl + AntrS - 4np] (2.9) 

The primes here denote derivatives with respect to x. The equation (2.3) is the Hamiltonian 
constraint while (2.5) is the momentum constraint. The constant mean curvature condition 
leads to the lapse equation (2.4). Equation (2.6) is a consequence of the choice of spatial 
coordinate condition. Equations (2.7) and (2.8) come from the definition of the second 
fundamental form and (2.9) is the one independent Einstein evolution equation which 
exists in this situation. To give the definition of the matter quantities occurring in (2.9) 
it is convenient to introduce a (locally defined) orthonormal frame. Let cq be the future- 
pointing unit normal to the hypersurfaces of constant t. Let ei be a unit vector tangent to 
these hypersurfaces which is normal to the surfaces of constant t and x. Complete eo and ei 
to an orthonormal frame by adding vectors 62 and 63. Then p = Tq^CoCq, j = —TctpeQe^, 
Sl = TapCiCi and trS = Tap{efei + 6263 + 63 Cg). 

Malec and O Murchadha [11] have written the constraints in an alternative form 
which turns out to be very useful in certain circumstances. (In fact they only consider 
the spherically symmetric case but the extension to plane and hyperbolic symmetry is 
straightforward.) They use as fundamental variables the expansions 9 and 9' of the families 
of null geodesies orthogonal to the orbits. (The prime in 9' is not a derivative.) In terms 
of the above coordinates these are given explicitly by 

9 = 2A-^A' + K -t 

(2.10) 
9' = 2A-^A' -K + t 

The sign conventions for 9 and 9' are those of [11] but the sign convention for the second 
fundamental form is the opposite of that used in [11]. The normalization of 9 and 9' 
depends on a choice of Cauchy hypersurface. Equation (2.10) applies to the case where 
the hypersurface chosen is a level surface of the coordinate t. However the analysis of [11] 



applies to any Cauchy hypersurface 5" compatible with the product decompostion of M 
in the sense that it is a union of surfaces of symmetry. The mean curvature of a Cauchy 
hypersurface of this kind will be denoted by H and is in general not constant. Let / be a 
proper length parameter along a curve in S orthogonal to the surfaces of symmetry. When 
expressed in terms of 9, 6' and / the constraints (2.3) and (2.5) become: 

die = -8n{p - j) - 16*2 -9H + e{aA)-^ 
diO' = -8n{p + j)-^e'^ + e'H + e{aA)-^ 

For the following discussion it is useful to introduce the area radius r = aA. The idea of 
[11] is to use the equations (2.11) to obtain bounds for the quantities r9 and r9'. It is 
assumed that the dominant energy condition holds so that \j\ < p. From (2.11) 



di(re) = -8nr(p - j) - — (^V^ - 4e + AOHr'^ + OriOr - O'r)) 

Ar 

dArO') = -8nr(p + j) - —(9'^r^ - 4e - 49' Hr'^ + e'riO'r - Or)) 

Ar 



(2.12) 



Consider now one particular Cauchy hypersurface 5" which is a union of surfaces of sym- 
metry. Denote the maximum value attained by r6 and r6' on this hypersurface by M_|_ and 
the minimum by M_. Let xq be a point where M_|_ is attained and suppose without loss 
of generality that 6{xq) > 9'{xo)- Since xq is a critical point of r9 it follows from (2.12) 
that at that point 

e^r^ + {AHr) (Or) - 4e < (2.13) 

Working out the roots of the corresponding quadratic equation then shows that 

-2{\Hr\ + \/e + H^r^) < M+ < 2{\Hr\ + y^e + H^r^) (2.14) 

The same inequality holds for M_. For each symmetry type this implies that r9 and r9' 
can be bounded in terms of H and r. An important quantity in the following is the mass 
function m which is defined by 

e-2m/r= ^r^ee' (2.15) 

The boundedness result just obtained shows that the following holds: 
Lemma 2.2 Let {M,g) be a surface symmetric spacetime which is foliated by compact 
CMC hypersurfaces with the mean curvature varying in a finite interval (^1,^2). If the 
dominant energy condition holds and r is bounded then 2m/r is bounded. 

There is another formulation of the equations which is also useful. This is based on the 
fact that the field equations can be written as equations on the two dimensional manifold 
of symmetry orbits. In the following lower case Latin indices take the values and 1 and 
are used to express tensor equations on this quotient manifold. In particular gab and Tab 
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denote the tensors on the quotient manifold naturaUy related to the spacetime metric and 
the energy-momentum tensor. The equations of importance in the following are 



m 



'^a'^br = —gab - 4.7ir{Tab " tvTgab 
Vafn = 4nr''{Tab - trTgab)V\ 



(2.16) 
(2.17) 



Note that these equations do not contain e explicitly. The expression for m in this formu- 
lation is 

m = ir(e-VWar) (2.18) 

The following is a generalization of an argument of Burnett[2]. 

Lemma 2.3 Let (M, g) be a surface symmetric spacetime and let S" be a compact Cauchy 

hypersurface of the form S" x F for a curve 5" in R x S"^. Let rmin and rmax denote the 

minimum and maximum values of r on S" respectively. Define R^ to be rmin, or — rmax 

for e equal to 1, or —1 respectively. Then 2m > R^ on S. 

Proof Let U be the set of x & S where 2m < R^. Then V"r is spacelike on U. Let Sa be 

the projection of the gradient of r onto 5". Then 



s'^Vafn = {Tab - trTgab)s''V\ 



(2.19) 



If the dominant energy condition holds then the right hand side of (2.19) is non-negative. 
Thus on U the mass m increases in the direction in which r increases. Now the restriction 
of r to S" cannot have a stationary point in U. This means in particular that S\U is 
non-empty. If xq is a point of U then moving from xq in the direction of increasing r leads 
to an increase in m. Eventually a point of the boundary of U must be reached. At that 
point m = Re/2. Thus 2m{xo) < R^. Moving away from xq in the opposite direction gives 
the reverse inequality. Hence 2m = R^ on U. However this contradicts the definition of U 
unless U is empty. 

In the spherical case Lemma 2.3 implies that the minimum of 2m, on S is greater 
than or equal to the minimum of r there. For plane symmetry it implies that the mass 
is no n- negative. The latter conclusion can be strengthened using (2.11) to get a kind of 
positive mass theorem. 

Lemma 2.4 Let (M, g) be a plane symmetric spacetime which satisfies the dominant 
energy condition. Then if to = at any point the spacetime is flat. 

Proof Consider any compact Cauchy hypersurface S of the form S x F. If to vanishes at 
some point xq of S then (for e = 0) either 9 or 9' must vanish there. Suppose without loss 
of generality that it is 9. Equation (2.11) gives 



di9 = -{H + 19)9 - 8tv{p - j] 



(2.20) 



This equation is similar to one which arises in a similar context for Gowdy spacetimes[5] 
and can be treated in exactly the same way. In fact if / is an arc length parameter which 
is zero at xq then the solution of (2.20) is 



^(/) = -87r / {p-j){u)e^p 



h-H-\9){v)dv 

J u 



du 



(2.21) 



From this formula it is clear that 6 is everywhere non-positive and that it can only become 
zero for some positive / if p — j vanishes identically on the interval [0,/]. In that case 9 
also vanishes on that interval. Since 9{l) is a periodic function it follows that it must be 
identically zero and that p = j everywhere on S. The mass is also zero on S. When 9 is 
zero the rate of change of r along S is given by 9' . Since the restriction of r to S must 
have a critical point somewhere, 9' must vanish somewhere. Applying to 9' the argument 
previously applied to 9 shows that 9' = and that p = j = 0. By the dominant energy 
condition p = implies the vanishing of the whole energy-momentum tensor on S. When 
the dominant energy condition holds the vanishing of the energy-momentum tensor on a 
Cauchy hypersurface implies that it vanishes everywhere. Thus the spacetime is vacuum. 
This in turn implies that m is zero on the whole spacetime. It follows that 9 and 9' are 
identically zero and that r is constant. In vacuum the Gaussian curvature of the metric 
Qab is K = r~^Ar and so in the present case K = and gab is flat. It is easily seen that 
this and the constancy of r imply that the spacetime is flat. 

In the case of hyperbolic symmetry the following analogue of Lemma 2.4 holds. 
Lemma 2.5 Let {M,g) be a spacetime with hyperbolic symmetry which satisfles the 
dominant energy condition. Then V"r is timelike. 

Proof V"r must be timelike somewhere on a symmetric Cauchy surface, as shown in the 
proof of Lemma 2.3. Hence it suffices to show that VafV^r never vanishes. This follows 
immediately from the analogue of (2.21). 

Non-flat spacetimes with plane symmetry and all spacetimes with hyperbolic symme- 
try have the property that the gradient of r is either everywhere past-pointing timelike or 
everywhere future-pointing timelike. Spacetimes where the first possibility is realized may 
be called 'expanding' models those where the second is realized 'contracting' models. This 
terminology is justified by the following considerations. If the gradient of r is past pointing 
then 9 is positive and 9' is negative. From (2.10) then follows that t < K everywhere. If t 
were non- negative then this would mean that \K\ was everywhere greater than |t|. However 
this is inconsistent with the existence of a compact CMC hypersurface, as can be seen by 
integrating the Hamiltonian constraint (2.3). Thus if the gradient of r is past-pointing in 
a region foliated by compact CMC hypersurfaces then t must be negative. Similarly, if the 
gradient of r is future-pointing, t must be positive. By possibly replacing t by —t it can 
be assumed without loss of generality that the model is expanding i.e. that the r increases 
monotonically with t along any causal curve. In that case t < in the whole spacetime so 
that ti < 0. In this case it will be assumed that in fact t2 < 0. 

The information obtained so far in this section implies bounds for various geomet- 
rical quantities in a surface symmetric spacetime without centre defined on a finite time 
interval (^1,^2). Now as many other quantities as possible will be bounded. Suppose that 
the energy-momentum tensor satisfies the non-negative pressures condition. Then (2.17) 
implies that when Va^ is timelike the rates of change of r and m along an integral curve 
of VaT" have opposite signs. It follows that in the cases e = and e = — 1 the radius r is 
bounded below by a positive constant and the mass bounded above on the interval (^3, ^2) 
while the radius is bounded above and the mass bounded below on the interval (^1,^3), 
where ^3 is any time with ti < ts < ^2. If e = or if e = — 1 and it is assumed that m is 



positive for t = t^ then the mass is bounded below by a positive constant on (^1,^3). In 
the case e = 1 results of Burnett [2] show that the radius is bounded above and the mass 
bounded below by a positive constant on both these intervals. Using the upper bound for 
2m/ r obtained earlier it can be seen that on an interval where the radius is bounded above 
and the mass bounded below by a positive constant the mass is bounded above and the 
radius is bounded away from zero. 

An estimate for the lapse function a can be obtained from (2.4). Considering a point 
where a attains its maximum and using the fact that p + trS > shows that a < 3/t^. 
Hence if ti 7^ and ^2 7^ it follows that a is bounded on (ti, 12). An interval which satisfies 
ti 7^ and t2 7^ and where r is bounded above and m is bounded below by a positive 
constant will be called admissible. The volume of the slice t = const, is V{t) = C J„ a'^A'^ 
and its time evolution is given by 

dV/dt = -Ct / aa^A^ (2.22) 

Jo 

In the case of an admissible interval this shows that V{t) and its inverse are bounded. Now 
V = a~^ Jq r^ and so if ti 7^ 0, t2 7^ there are positive constants Ci and C2 such that 
Ci ^ a < C2 and Ci < A < C2. Now integrate the equation (2.3) from to 27r. There 
results the inequality 

r 16npA'>/^ < 2 r A-^/^a-^ + | /'" t^A^/^ (2.23) 

Jo Jo Jo 

which implies that on an admissible interval j^ p is bounded. The dominant energy 

condition then gives bounds for j j and /trS. The bounds for r6 and r9' now show that 

A' and K are bounded on any admissible time interval. Integrating (2.6) over the circle 

allows at to be bounded and then (2.6) itself gives a bound for 13' . A bound for [3 can be 

deduced using the condition J/3 = 0. Equation (2.8) gives a bound for At. Integrating 

equation (2.4) from a point where a' = and using (2.23) provides a bound for a' . 

Theorem 2.1 Let a solution of the Einstein equations with surface symmetry be given 

and suppose that when coordinates are chosen which cast the metric into the form (2.1) 

with constant mean curvature time slices the time coordinate takes all values in the finite 

interval (ti,t2). Suppose further that: 

i) the dominant energy and non-negative pressures conditions hold 

ii) neither ti or ^2 is zero 

iii) if e is or 1 then ti < 

iv) if e = — 1 then the mass function is positive on the initial hypersurface. 

Let ta satisfy ti < ^3 < ^2- Then the following quantities are bounded on (^1,^3): 

a, a', A, A-\ A', K, (3, a, a'^, dta (2.24) 

dtAK',/3' (2.25) 



3. The matter fields 

Theorem 2.1 provides some information on the boundedness of certain geometrical 
quantities in a spacetime with spherical, plane or hyperbolic symmetry without any as- 
sumptions on the matter content except the dominant energy and non-negative pressures 
conditions. To get further bounds and hence to proceed towards showing that the space- 
time can be extended it is necessary to use the matter field equations. In the following 
two examples will be treated, namely the collisionless gas and the massless scalar field. 

The collisionless gas is described by a distribution function / which is a non-negative 
real-valued function on the mass shell. It is supposed to satisfy the Vlasov equation which 
in the class of spacetimes considered here takes the form: 



df 



dt 



^+ aA-'^-p]^ + 



df 



— a 



A- 



dx 

v^ 1 

V- 2 



^';;o + o(^-^) 



-A-^av^ + aKv^ + aA-'^A 
Bdf 



3^2 



^{v'y + iv'') 






(3.1) 



dv 



B 







Here the mass shell has been coordinatized using components in an orthonormal frame, 
where the first vector in the spatial frame is proportional to d/dx. The component v^ 
is then given by the expression a/1 + (v^)'^ + (f^)^ -|- (v^)'^. The upper case Latin indices 
take the values 2 and 3. The distribution function depends on t, x, v^, v"^ and v^ . In fact 
the symmetry requires that its dependence on the last two quantities is only a dependence 
on the combination {v"^)"^ + (t'"^)^. This will be assumed for the initial data and is then also 
satisfied by the solution. The initial data is assumed to be compactly supported and then 
the solution has compact support at each fixed time. That the initial value problem for 
the Vlasov-Einstein system is well posed was shown by Choquet-Bruhat[3]. The matter 
quantities occurring in the field equations are given by 



P 

J 

Si 

trS 



fv^dv 
fv^dv 

fiv^f/v'^dv 
/[(z;°)2 - V^/v'^dv 



(3.2) 



For a solution which evolves from initial data given at t = to let 



P{t) = supjlfl : /(s, X, f ) 7^ for some (s, x, v) with s G /} 



(3.3) 



where / is the interval [to,t] if t > to and the interval [t, to] if t < to. The maximum of 
/ is time independent and so all the matter quantites defined in (3.2) can be bounded by 
C(l + P(t))'^. The quantity P(t) itself can be controlled by studying the characteristics of 
the equation (3.1) since this equation says that / is constant along these characteristics. 
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It follows that P(t) will be bounded on a given interval provided the coefficients in (3.1) 
are bounded. The geometrical quantities which occur in these coefficients are a, A~^, /3, 
K and t. Theorem 2.1 shows that all of these are bounded on an admissible time interval. 
Considering a point where a attains its minimum on a given hypersurface t=const. leads 
to a bound for a~^. Hence the following is obtained. 

Theorem 3.1 Let a solution of the Vlasov- Einstein system with surface symmetry be 
given and suppose that when coordinates are chosen which cast the metric into the form 
(2.1) with constant mean curvature time slices the time coordinate takes all values in the 
finite interval (ti, ^2)- Suppose further that conditions (ii)-(iv) of Theorem 2.1 are satisfied. 
Then all the quantities in (2.24) and (2.25) are bounded (^1,^3), as are P, a~^ and 

pJ.SltrS (3.4) 

Notice that while Theorem 2.1 is to a large extent independent of the matter model 
used this is not the case for Theorem 3.1. It is probable that the analogous statement would 
be false if the coUisionless matter was replaced by dust since shell-crossing singularities 
would presumably provide counterexamples. It will now be shown by induction that under 
the hypotheses of Theorem 3.1 all derivatives of the solution are bounded. 
Lemma 3.1 If the hypotheses of Theorem 2.1 are satisfied and if all derivatives with 
respect to x of order up to n of the quantities in (2.24) and (3.4) are bounded then all 
derivatives with respect to x of order up to n + 1 of the quantities in (2.24) are bounded. 
Proof In the following D^ denotes a derivative with respect to x. Note first that the 
boundedness of the derivatives with respect to x up to order n + 1 of the quantities a, 
A, A~^, a and a~^ follows immediately from the hypotheses of the lemma. Equations 
(2.3)-(2.6) can be solved for the quantities A" , a" , K' and (3' . Differentiating the resulting 
equations n times with respect to x allows D2~^^{A'), D2~^^{a'), D^'^^K and D'^'^^jS to be 
bounded. 

Lemma 3.2 If the hypotheses of Lemma 3.1 are satisfied by a solution of the Vlasov- 
Einstein system then all derivatives with respect to x of order up to n + 1 of the quantities 
in (3.4) are bounded. Moreover all derivatives of / of order up to n + 1 with respect to x 
and V are bounded. 

Proof The hypotheses imply that the coefficients of the characteristic system are bounded 
together with their derivatives with respect to x up to order n + 1. Differentiating this 
system n + 1 times with respect to x and v gives an inhomogeneous linear system of 
ordinary differential equations for the derivatives of order n + 1 of the unknowns. The 
coefficients of this system are bounded, as long as attention is confined to the support of 
/. Hence these derivatives are bounded on the support of /. It follows from this that the 
derivatives of / up to order n + 1 are bounded. Differentiating (3.2) then gives the desired 
conclusion for the quantities in (3.4). 

Lemma 3.3 If the hypotheses of Theorem 2.1 are satisfied by a solution of the Vlasov- 
Einstein system and if all derivatives of the quantities in (2.24) and (3.4) of the form 
D^D'^ with n arbitrary and k < m are bounded and if all higher derivatives of / with 
at most m time derivatives are bounded then the derivatives of the form -D[" -D" of the 
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quantities in (3.4) are bounded. Moreover the higher derivatives of / with at most m + 1 
time derivatives are bounded. 

Proof Use the Vlasov equation to bound an extra time derivative of / and substitute the 
result into (3.2). 

Lemma 3.4 If the hypotheses of Lemma 3.1 are satisfied and if: 

(i) aU derivatives of the quantities in (2.24) and (3.4) of the form D^D"^ with n arbitrary 

and k < m are bounded 

(ii) aU derivatives of the quantities in (3.4) of the form -D["^ D"^ are bounded 

then the derivatives of the quantities in (2.24) of the form D^^^D"^ are bounded. 

Proof The conclusion for a and a~^ follows immediately from the assumptions. The 

conclusion for A, A~^, A' and K follows from equations (2.8) and (2.9). Now differentiate 

(2.4) n times with respect to x and m + 1 times with respect to t. The result is an equation 

of the form 

(D'^+^D^a)" + A-^A'iD'^+^D^a)' 

= A^iliK - it)2 + 1^2 ^ 4^(^ ^ trS)]L>r^^L>> + B 

where the remainder term B is bounded. Examining the points where D^^ D'^a has 
maximum modulus gives \DY^^^D'^a\ < 3A~'^B/t^. Next, the equation (2.7), integrated 
with respect to x, allows the conclusion to be obtained for dta. Finally, the conclusion for 
j3 follows from (2.6). 

Putting together the conclusions of Theorem 3.1 and Lemma 3.1-3.4 we see that under the 
hypotheses of Theorem 3.1 all derivatives of all metric coefficients and of the distribution 
function are bounded. 

When all the derivatives of the solution are bounded on a given interval, it can be 
extended smoothly to the closure of that interval. There results a new initial data set 
and applying the local existence and uniqueness result discussed in Appendix 1 gives an 
extension of the solution to an interval which strictly contains the original one. Thus we 
have the following theorem. 

Theorem 3.2 Let (M, g, f) be a C°° solution of the Vlasov-Einstein system with surface 
symmetry which is the maximal globally hyperbolic development of data given on a hy- 
persurface of constant mean curvature Hq. Then: 

1. If e = 1 and Hq = then the whole spacetime can be covered by a foliation of CMC 
hypersurfaces where the mean curvature takes all real values. 

2. Ife = lore = and Hq < then the part of the spacetime to the past of the initial 
hypersurface can be covered by a foliation of CMC hypersurfaces where the mean curva- 
ture takes all values in the interval (— cxd, Hq] 

3. If e = — 1, Hq < and the mass function is positive on the initial hypersurface then the 
part of the spacetime to the past of the initial hypersurface can be covered by a foliation 
of CMC hypersurfaces where the mean curvature takes all values in the interval (— oo, Hq] 

Now another example, the massless scalar field, will be discussed. The massless 
scalar field is described by a real-valued function (f) satisfying Vq,V"(/) = 0. The energy- 
momentum tensor is 

Tc.(3 = Va<^V/3(/) - I {\/^(PV^(P)gc.0 (3.5) 
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The dominant energy condition is satisfied but the non-negative pressures condition does 
not hold for a general (j). Hence Theorem 2.1 does not apply to this case. However, as 
has been remarked by Burnett [2], it is true that T^px^^x^ > for spacelike vectors x" 
orthogonal to the surfaces of symmetry and it turns out that this fact and the condition 
that p + trS > are the only ones which are needed in the proof of the theorem. They are 
satisfied by the scalar field (see below). There is also a potential problem with applying 
the analysis of Appendix 1 in this case, since there the non-negative pressures condition 
was also used. It was needed for the argument using the implicit function theorem but 
only in the case that the initial hypersurface is maximal. Hence the analogue of Theorem 
2.1 holds for the scalar field under the extra hypothesis that the interval (^1,^2) does not 
contain zero. 

To proceed further it is useful to introduce the null vectors 



dt dxj dx 

dt dx J dx 



(3.6) 



Let (j)+ = e+(j) cind (p- = e-4>. Then the wave equation for (p can be written as 

ae+{(f)-) = at{(j)+ + </)_) + {A'a + 2aA~^A'){(f)+ - (/>_) + a[e_, e+](/) (3.7) 

ae_((/)+) = at{(p+ + (p-) + {A'a + 2aA~^A'){(p+ - (p_) + a[e+, e_](/) (3.8) 

Now Q;[e_|_,e_] = 6+e+ + 6_e_, where the coefficients 6_|_ and 6_ are polynomials in the 
quantities (2.24) and (2.25). In addition, the definitions of e-|- and e_ imply that 

d(l)/dt-l3d(l)/dx = ^a{e+ + e-) (3.9) 

The matter quantities occurring in the field equations are given by 

p=l{<Pl + <pl) 

Sl = l{cPl + cpl) (3.10) 

trS = i((/>+ + (^_)2 + i<^+</,_ 
p + trS = i((/)+ + (^_)2 

Let 

$(t) = \m)\\o. + u+moo + u-moo (s.n) 

The equations (3.7)-(3.9) and the boundedness of the quantities (2.24) and (2.25) imply 
that 

$(t) < $(0) + C f ^{s)ds (3.12) 

Jo 
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Hence (/>, (/>_(_ and (j)- are bounded on an admissible interval. It then follows from (3.10) that 
the matter quantities in the field equations are bounded. Hence the analogue of Theorem 

3.1 with the Vlasov-Einstein system replaced by the Einstein-scalar system and P replaced 
by $ holds provided the interval (^1,^2) does not contain zero. 

Lemma 3.5 If the hypotheses of Lemma 3.1 are satisfied by a solution of the Einstein- 
scalar system then all derivatives with respect to x of order up to n + 1 of the quantities 
in (3.4) are bounded. Moreover all derivatives of (j) of order up to n -|- 1 and of dt4> up to 
order n with respect to x are bounded. 

Proof Note first that the results of Lemma 3.1 can be improved slightly by using (2.6) 
and (2.8) to bound the spatial derivatives of order up to order n + 1 of the quantities 
13' and dtA. Next, differentiating the equations (3.7)-(3.9) gives an inhomogeneous linear 
hyperbolic system for L>"+^(/)_|_, D'^^^(p- and W^^^cp with a bounded right hand side. This 
gives the desired bounds for derivatives of (j). Putting these into (3.10) gives the bounds 
for the components of the energy-momentum tensor. 

Lemma 3.6 If the hypotheses of Theorem 2.1 are satisfied by a solution of the Einstein- 
scalar system and if all derivatives of the quantities in (2.24) and (3.4) of the form D^D^ 
with n arbitrary and k < m are bounded and if all higher derivatives of (p with at most 
m + 1 time derivatives are bounded then the derivatives of the form D^~^ D"^ of the 
quantities in (3.4) are bounded. Moreover the higher derivatives of (p with at most m + 2 
time derivatives are bounded. 
Proof This result follows immediately from the equations (3.7)-(3.9) and (3.10) 

Putting together the analogue of Theorem 3.1 for the scalar field and the Lemmas 
3.1, 3.4, 3.5 and 3.6, we see that under the hypotheses of the analogue of Theorem 3.1 and 
assuming that Hq 7^ all derivatives of all metric coefficients and of the scalar field are 
bounded. If Hq = then we still get a neighbourhood of the initial hypersurface foliated 
by CMC hyper surfaces. If there exists a hypersurface belonging to this foliation with non- 
zero mean curvature both to the past and to the future of the initial hypersurface then 
the problem of getting bounds is reduced to the case Hq 7^ 0. If all the CMC hypersur- 
faces belonging to the foliation which are to the past, say, of the initial hypersurface are 
maximal then the second fundamental form and the Ricci tensor contracted twice with the 
normal vector are both zero on that region. The latter implies that the gradient of (p is 
tangent to the CMC hypersurfaces. Under these circumstances the wave equation reduces 
to the Laplace equation and (p must be spatially constant. If (p varied from one spacelike 
hypersurface to the next then V°'(p would be timelike, contradicting what has been said al- 
ready. Hence (p is constant, the energy-momentum tensor is zero and the vacuum Einstein 
equations are satisfied. These considerations show that the following analogue of Theorem 

3.2 holds for the scalar field: 

Theorem 3.3 Let {M,g,(p) be a C°° solution with surface symmetry of the Einstein 
equations coupled to a massless scalar field which is the maximal globally hyperbolic 
development of data given on a hypersurface of constant mean curvature Hq. Then: 

1. If e = 1 and Hq = then the whole spacetime can be covered by a foliation of CMC 
hypersurfaces where the mean curvature takes all real values. 

2. Ife = lore = and Hq < then the part of the spacetime to the past of the initial 
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hypersurface can be covered by a foliation of CMC hypersurfaces where the mean curvature 
takes all values in the interval (—00, Hq] 

3. If e = —1, Hq < and the mass function is positive on the initial hypersurface then the 
part of the spacetime to the past of the initial hypersurface can be covered by a foliation 
of CMC hypersurfaces where the mean curvature takes all values in the interval (—00, Hq] 

Appendix 1 

The purpose of this appendix is to prove Lemma 2.1 and a local existence and unique- 
ness theorem for the equations (2.3)-(2.9). Let (M, got/3) be a surface symmetric spacetime, 
as defined in Section 2. Thus, in particular M is of the form Hx S^ x F. Let M be the uni- 
versal cover of M and let ^a/3 the pull-back of g^p to M. Suppose that (M, gotp) contains a 
CMC Cauchy hypersurface 5". It follows from the fact [10] that all globally defined Killing 
vectors must be tangent to a compact CMC hypersurface that a compact CMC Cauchy 
hypersurface 5" in a spacetime with spherical or plane symmetry is of the form S x F. For 
in that case there are enough global Killing vectors to generate the surfaces of symmetry. 
This is not true in the case of hyperbolic symmetry. It is a standard fact that a neigh- 
bourhood of a compact CMC hypersurface can be foliated by compact CMC hypersurfaces 
unless the original hypersurface is such that its second fundamental form vanishes and 
the Ricci tensor contracted twice with the normal vector is zero. Furthermore the mean 
curvature of these hypersurfaces can be used as a time coordinate in this neighbourhood. 
Even if this condition fails there is still a neighbourhood of the initial hypersurface foliated 
by CMC hypersurfaces, although in that case the mean curvature cannot be used as a time 
coordinate [1]. When the dominant energy and non-negative pressures conditions are satis- 
fied the condition can only fail if the spacetime is vacuum on S. From what was said earlier 
the CMC hypersurfaces must be symmetric in the cases of spherical and plane symmetry. 
In fact they must also be symmetric in the hyperbolic case. To see this note that the 
existence of these hypersurfaces is proved by using the inverse function theorem. However 
it is possible to apply the inverse function theorem in the class of symmetric deformations 
of the initial hypersurface (provided this initial hypersurface is itself symmetric) and then 
the CMC hypersurfaces obtained are by construction symmetric. Consider now the excep- 
tional case where the spacetime is vacuum and the second fundamental form of the initial 
hypersurface vanishes. In the cases of spherical and hyperbolic symmetry integrating (2.3) 
from to 27r gives a contradiction. In the plane symmetric case the vanishing of the second 
fundamental form implies that spacetime is fiat. Hence unless the spacetime is fiat the 
mean curvature can be used as a time coordinate t in a neighbourhood of S. The inverse 
image of t under the projection p : M -^ M will also be denoted by t. 

It is elementary to see that the metric gab of the hypersurface t = const, in M can be 
written locally in the form 

A^dx^ + B^dE'^ (Al.l) 

where (iE^ is a metric of constant curvature. What is less clear is that that this can be done 
globally in such a way that x G [0, 27r] and the functions A and B are 27r-periodic. Consider 
one of the hypersurfaces t = const, in M. Let 7 be a geodesic in this hypersurface which 
starts orthogonal to one of the group orbits Oi at a point p. It continues to be orthogonal 
to the orbits. After a finite time it must hit an orbit O2 which projects to the same orbit 
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in M a,s Oi. Suppose it meets O2 at a point q. Let q' be the unique point of Oi which 
projects to the same point of M as q. Any isometry which fixes p must fix q' . It follows 
in the plane and hyperbolic cases that p = q' and in the spherical case that either p = q' 
or p and q' are antipodal points on the sphere. If p = q' then A and B can be made 
27r-periodic, as desired. In the case where the antipodal map occurs the same thing can 
be arranged by allowing x to go twice around the circle. 
Let 

r.271- 



B{x)/A{x)dx 



(A1.2) 



a — 27T 
Then the new coordinate x' defined by 

x' = a I B{x)/A{x)dx (Al.3) 



satisfies x(0) = and x{2ti) = 2tt. Define A' by the relation A'{x')^{x) = A{x). After 
transforming to the new coordinate and dropping the primes the metric takes the form 

A'^{dx^+a^dJ:^) (AlA) 

where A is a positive function of x with A{0) = A{27r) and a is a constant. Doing this 
construction on each hypersurface of constant time gives the coordinate system whose 
existence is asserted by Lemma 2.1. 

Consider now the question of local existence and uniqueness of solutions of equations 
(2.3)-(2.9) with given initial data on a hypersurface t = const.. In order to have a well-posed 
initial value problem it is necessary to have some matter equations such that the resulting 
Einstein-matter system has a well-posed Cauchy problem in the context of C°° data and 
solutions. More precisely we assume that the solution of the initial value problem for the 
reduced equations in harmonic coordinates exists and is unique so that the general theory 
of the maximal Cauchy development [4] can be applied. (If it were desired to consider 
gauge theories, where solutions are only unique up to gauge transformations, then some 
more work would be required.) It is not obvious that this theory applies to kinetic theory 
models, where the matter fields are defined on the mass shell rather than on spacetime. 
However the analogous results do hold in that case [3]. An initial data set consists of 
periodic functions A and K, a constant a and matter data which satisfy the constraint 
equations (2.3) and (2.5). The matter data are assumed to have the necessary symmetry 
properties. These properties are most easily expressed on the covering manifold. The 
data set on the covering manifold has a maximal Cauchy development on the manifold M. 
The maximal Cauchy development inherits the symmetries of the data and so the original 
initial data set has a surface symmetric Cauchy development. In this surface symmetric 
spacetime coordinates can be introduced as above. Thus a solution of equations (2.3)-(2.9) 
(and the matter equations) on some interval (ti,t2) is obtained. It remains to show that 
solutions of these equations are uniquely determined by initial data. Suppose there exist 
two solutions with the same initial data on the interval (ti, ^2). Then by the general theory 
of the Cauchy problem there must exist embeddings (pi and (p2 of M into the maximal 
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Cauchy development of the given initial data set such that (pi is a matter preserving 
isometry for the first solution and (/>2 a matter preserving isometry for the second. The 
uniqueness of compact constant mean curvature hypersurfaces implies that the images of 
any hypersurface of constant t under (pi and (p2 are identical. In particular this means 
that the images of M under (pi and (p2 are identical, so that there exists a diffeomorphism 
(pi2 : M ^ M such that (pi = (pi2 o(p2- The diffeomorphism (pi2 maps the one solution into 
the other and preserves the hypersurfaces of constant time. Suppose temporarily that the 
initial data do not have Robertson- Walker symmetry. Then a unique two-plane is defined 
by the isotropy group of the universal cover. This plane must be preserved by (/)i2 as must 
its orthogonal complement. Thus, if (pi2 is written in terms of coordinates adapted to 
the first metric in the form {t,x,y) i-^ {t',x',y'). Then t' = t and y' depends only on t 
and y. By composing with an isometry it can be reduced to the identity on the initial 
hypersurface. The form of the shift vector then implies that it is the identity everywhere. 
In a given spacetime the coordinate x is defined up to a translation. Hence x' — x + c and 
since (pi2 is the identity on the initial hypersurface it follows that the two solutions are 
identical. In the case where the data have Robertson- Walker symmetry the solutions must 
have Robertson- Walker symmetry and in that case uniqueness for the reduced equations 
is obvious. 

Appendix 2 

In this appendix the vacuum solutions with the symmetry properties considered in 
this paper will be determined. This is done using the following lemma: 
Lemma A2.1 Consider the ordinary differential equation d'^u/dx'^ = f{u), where / : 
(0, oo) ^ R is Lipschitz. Suppose that /(wq) = for some uq, f{u) < for < u < uq 
and f{u)>0 for u > uq. Then any periodic solution is constant. 

Proof Let w be a periodic solution. By periodicity there exists a point xq where d'^u/dx'^ 
vanishes. At that point u = uq. If du/dx{xQ) is positive then it is easy to show that 
du/dx remains positive for x > xq, contradicting periodicity. Similarly the assumption 
du/dx{xo) < leads to a contradiction. Hence in fact du/dx{xo) = 0. By uniqueness for 
solutions of the ordinary differential equation it follows that u is constant. 

Consider now vacuum solutions of equations (2.3)-(2.9). The momentum constraint 
can be solved explicitly, giving K — ^t = CA~^ for some constant C. Substituting this 
into the Hamiltonian constraint gives: 

(AV2)" = -3_c^A-y^ + ^t2^^/2 ^ iea-2A-i/2 (A2.1) 

It can be checked straightforwardly that this ordinary differential equation for A^''^ satifies 
the hypotheses of the lemma and so A is constant. Then the same lemma may be applied 
to the lapse equation to show that a is constant. The constancy of A implies that of K 
and the equation for P then gives P = 0. Hence every vacuum solution of equations (2.3)- 
(2.9) is spatially homogeneous. These solutions will now be identified with known exact 
solutions. This will be done by examining the Cauchy data on one spacelike hypersuface. 
Suppose that constants t, a and K are given and satisfy the following sign condition, which 
is necessary for the constraints to have a solution: 
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(i) if e = 1 then ^{K - ^t)^ - ft^ > 
(ii) if e = then f (K - ^tf - ft^ = 
(iii) if e = -1 then f (K - ^tf - ft^ < 

Suppose that t = 0. Then the sign condition is incompatible with e = —1. It is only 
compatible with e = if K = 0. In that case the data give rise to fiat space, identified in 
a simple way. If e = 1 then the Hamiltonian constraint can be solved for A in terms of a 
and K. For t y^ the sign condition can readily be analysed by dividing the expression of 
interest by t^ and studying the resulting quadratic expression in K/t. 

The case e = is the simplest. There are two possible values for K/t, namely —1/3 
and 1. These solutions of the constraints can be realized by the r=const. hypersurfaces 
in the Kasner solution 

-dr^ + b^T^Pdx^ + T^-P{dy^ + dz^) (A2.2) 

where p=— l/3orp = l and 6 is a positive constant. In the case e = 1 the quantity 
K/t takes all values in the intervals (— oo, — 1/3) and (l,oo) and these solutions of the 
constraints can be realized by the r=const. hypersurfaces in the following metric, which 
is obtained by identifying the part of the Schwarzschild solution inside the horizon: 

-(2m/r - 1)-Mr2 + }p{2m/T - l)dx'^ + r^rfE^ (A2.3) 

Here dT? is the standard metric on the sphere. Similarly, the solutions with e = — 1 produce 
all values of K/t in the interval (—1/3, 1) and these solutions of the constraints can be 
realized by the r=const. hypersurfaces in the following pseudo- Schwarzschild metric: 

-(2m/r + l)-^rfr2 + }p{2m/T + l)dx'^ + r'^dT? (A2.4) 

In this case dT? is a metric of constant negative curvature on a compact manifold obtained 
by identifying the hyperbolic plane by means of a discrete group of isometries. The general 
theorems proved in this paper imply in particular that for m > the initial singularity in 
this solution, which occurs at t = is a crushing singularity. It is worth remarking that 
for m < the initial singularity, which occurs at t = —2m is also crushing, even through 
the theorems do not apply. 
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